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Abstract. - The analysis of dissipation and dephasing in driven mesoscopic devices requires a 
distinction between two notions of quantum irreversibility. One ("Loschmidt echo") is related 
to "time reversal", while the other is related to "driving reversal". In the latter context we 
define the time of maximum return ("compensation time") which generalizes the notion of 
"echo time". Non-perturbative features manifest themselves in the energy spreading process. 
This is demonstrated for the prototype random-matrix Wigner model, where the compensation 
time and the system response exhibit a non-universal scaling behavior. 



Quantum Irreversibility (QI) is a subject of a recent intensive research activity due 
to its relevance to quantum computing [2] . It turns out that the analysis of dissipation and 
dephasing in driven mesoscopic devices requires a distinction between two different notions 
of "irreversibility". One is based on the "piston model" paradigm (PMP), while the other 
is based on the "ice cube in a cup of hot water" paradigm (ICP). The latter notion has 
been adopted by the recent literature PG], and it is related to the studies of dephasing due 
to the interaction with chaotic degrees of freedom The PMP on the other hand has 
direct relevance to our recent studies 4 of the (non-perturbative) response of driven chaotic 
mesoscopic systems. 

In the PMP case we say that a process is reversible if it is possible to "undo" it. Consider 
a gas inside a cylinder with a piston. Let us shift the piston inside. Due to the compression 
the gas is heated up. Can we undo the "heating" simply by shifting the piston outside, back 
to its original position? If the answer is "yes", as in the case of strictly adiabatic process, 
then the process is said to be reversible. In the ICP case we consider the melting process of 
an ice cube. Let us assume that after some time we reverse the velocities of all the molecules. 
If the external conditions are kept strictly the same, we expect the ice-cube to re-emerge out 
of the water. In practice the external conditions (fields) are not exactly the same, and as a 
result we have what looks like irreversibility. 

It is also essential to define precisely what is the measure for QI. The prevailing possibility 
is to take the survival probability as a measure pQ. Another possibility is to take the energy 
spreading as a measure. (The energy spreading is defined as the square root of the energy 
variance.) The latter definition goes well with the PMP, and it has a well defined classical 
limit. 
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Thus the ICP is related to "time reversal" experiments, while the PMP is related to 
"driving reversal" experiments. Later we are going to put the focus on PMP irreversibility. 
A thorough understanding of the one-period driving reversal scenario is both important by 
itself, and also constitutes a bridge towards a theory for the response to periodic driving 0|. It 
should be realized that irreversibility in the PMP sense implies irreversible diffusion in energy 
space. This diffusion leads to irreversible energy absorption ("dissipation") in periodically 
driven mesoscopic devices One purpose of this Letter is to sharpen the notion "irreversible 
diffusion" . 

The structure of this Letter is as follows: First we discuss the mathematical formulation 
of the ICP and the PMP. Then we focus on QI in the PMP sense, which has not been 
done in the past as far as we know. In this context we define the time of maximum return 
("compensation time"), which constitutes a generalization of the notion of "echo time". The 
specific model that we study is the well known Wigner model UJ. We shall explain that 
this is the simplest model that contains the main ingredients of a physical problem. We 
define the notions of perturbative/non-perturbative regimes and in particular we explain the 
adiabatic/sudden limits. The nature of QI in the different dynamical regimes is found to be 
distinct. The theoretical considerations are supported by the analysis of extensive numerical 
simulations. The new insights regarding the analysis of wavepacket dynamics in energy space 
[S] are summarized in the last paragraph. 

Consider a system whose evolution is governed by the chaotic Hamiltonian 7i = 7i{Q a , P a ;x), 
where (Q a , P a ) is a set of canonical coordinates, and x = x(t) with < t < T/2 is a time de- 
pendent parameter. The corresponding evolution operator will be denoted by U[x\. In the ICP 
case the parameter x represents the external fields. The physical conditions that correspond 
to the "forward" evolution are represented by x = XA(t), while the physical conditions that 
correspond to the "backward" evolution are the time-reversal of x = XB(t). If Xs{t) — xa(£) 
then we have complete reversibility, otherwise we can define a generalized fidelity amplitude as 
F(T) = (4 , o|^[^b] _1 C^[^a]| v I'o), where is the initial preparation. The generalized fidelity 
amplitude F(T) is in fact [3] the Feynman- Vernon influence functional, and \F(t)\ < 1 is the 
so called dephasing factor. For constant xa and xb we can write U[xa] — exp(— iTHa/^K) 
and f [s;b] = exp(— iT7iB/2?i), leading to the currently prevailing definition of the fidelity, 
also know as Loschmidt echo, as a measure for quantum irreversibility £Q. 

To be more specific we assume that the variations of the external fields are classically 
small, so that we can linearize the Hamiltonian as follows: Ti = Hq + xW. Then, for constant 
perturbation x = ±e, we can write 7Ya = 7io+eW and H.b = Ti-o—eW, where e is the amplitude 
of the perturbation. Note that the prevailing definition of the fidelity Vi,e(T) — |F(T)| 2 , also 
known as Loschmidt echo, takes the form 

V LE (T) = \(y \ e ~^(- H °+^e-^ n °+^\y )\ 2 (1) 

In complete analogy with the ICP, we can define V(T) = |(* |[/[:EB(rev)][/[2; A ]|\l/ )| 2 . 
Here xb(jsv) corresponds to the reversed driving process, eg pushing the piston back to its 
initial location. For constant xa and xb we can write U[xa] — exp(— i^Ha) and J7[xB(rev)] = 
exp(— i-jH.B), and for the linearized Hamiltonian we get 

V(T) = |(* |e- l *(^ eW )e-^ Ho+eW >|*o)| 2 (2) 

If we take the PMP literally, then the parameter e represents either a small displacement of a 
piston, or its velocity. The latter interpretation emerges if we write the Schrodinger equation 
in the adiabatic basis (See eg Eq. (55) of Ref. [S])- Similarly, in case of an Aharonov-Bohm 
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ring system, the parameter e represents either the magnetic flux through a ring, or else it 
may represent the electro- motive-force (measured in Volts). The latter, by Faraday law, is 
the time derivative of the former. As a matter of terminology, note that if e is interpreted as 
"velocity", then what we call below "standard perturbative regime" is in fact the "adiabatic 
regime" (small e) . Later we also discuss the notion of " sudden regime" (large e) . 

At first sight the two definitions of V(T) look very similar: it is as if the role of 7io and W 
is interchanged. Also one may (wrongly) think that there is some kind of symmetry between 
the "adiabatic limit" (e — * 0), and the opposing "sudden limit" (e — > oo). This superficial 
"algebraic" symmetry of the definitions is misleading. One should be aware that in typical 
circumstances TLq has an unbounded spectrum, unlike the perturbation W. Specifically we 
refer here to quantized chaotic systems. In the standard representation TLq = E is a diagonal 
matrix whose elements are the ordered energies {E n }. We limit our interest to a classically 
small (but quantum mechanically large) energy interval where the mean level spacing A is 
roughly constant. Assuming that the classical motion is characterized by a finite correlation 
time Tci , it follows that the matrix W = B in the standard representation is a banded matrix 
which looks random |7j. Its bandwidth is A = 2Trh/T c i- The simplest model that captures this 
generic feature is the Wigner banded random matrix model (BRM) The non- vanishing 
couplings are within the band < \n — m\ < b. These coupling elements are zero on the 
average, and they are characterized by the variance a = ((IB^I 2 )) 1 / 2 . For sake of later 
convenience we define the bandwidth in energy units as = bA, and describe the band 
profile using a spectral function C(iS) that equals 2irha 2 / A for \Tioj\ < Af,, and zero otherwise. 
This is in fact the Fourier transform of C(r) = 26tT 2 sin(r/r c ;)/(r/r c /). The function C(r) 
describes the time autocorrelation of W in the Heisenberg picture. 

We are interested in the time evolution of the energy distribution. The corresponding 
time-dependent Schrodinger equation has been integrated numerically using the self-expanding 
algorithm of |Sj to eliminate finite-size effects (the remarkable accuracy is demonstrated in 
Fig. la). In all the simulations the initial state \*f?o) is assumed to be an eigenstate |m) of 
Ho- In the first half period (t < T/2) we have "conventional" wavepacket dynamics that we 
had studied in [Sj. What we want to discuss in this Letter is the dynamics for t > T/2, after 
the driving is reversed. The probability distribution after time t is Pt(n\m) — \(n\U(t)\m}\ 2 . 
Averaging over realizations of the Hamiltonian one obtains the average profile Pt(n — to). We 
characterize the evolving distribution using two different measures: The survival probability 
V(t) = V(t;T) = P t (r=0), and the energy spreading 5E(t) = 5E(t;T) = A x y/J2 r r2p t( r )- 
Thus the energy distribution at the end of the driving cycle is characterized by V(T) = V(T; T) 
and SE(T) = SE(T;T). 

Looking on Fig. lb and Fig. lc we see that it is natural to define a time t r of maximal 
return. This is the time, after the driving reversal, when we observe maximum compensation. 
The "compensation time" (t r ) is determined either by the minimum of the spreading 5E(t) 
or by the minimum of the total transition probability p(t) = (1 — V(t)). It is clear, by simple 
inspection of the plots in Fig. 1, that the time of "maximal return" is in general not t r = T 
but rather {T/2) < t r < T. Therefore we cannot regard this compensation as an "echo" 
phenomena (the term "echo" implies maximal return at t r — T). The dependence of t r on T 
for various values of e is presented in Fig. 2a and Fig. 2b. 

We want to understand the effect of driving reversal on the spreading/decay processes, 
and in particular to determine whether (or to what degree) the spreading/decay processes are 
reversible. The systematic strategy for the analysis of the dynamics consists of two stages: The 
first stage is to define the different e regimes in the corresponding parametric problem. As far 
as this Letter is concerned the important distinction is between the standard perturbative 
regime (e < A/a), the Wigner Lorentzian regime (A/a < e < 6 1 / 2 A/cr), and the non- 
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perturbative semicircle regime (e > 6 1 / 2 A/cr). The second stage is to consider the time 
dependent scenario separately in each regime. The driving is characterized by two parameters: 
One is the amplitude e, and the other is the period T. Thus we have a two dimensional 
parameter space (T, e) to study, in analogy with having (CI, A) space in case of periodic 
driving 0]. If the perturbation parameter e has the meaning of "velocity" , then the "standard 
perturbative regime" is in fact the " adiabatic regime" . By association one naturally wonders 
whether there is also a " sudden regime" . We shall come to this issue later on. 

It is natural to inquire whether perturbation theory is applicable for the analysis of the 
driving reversal scenario. Using first order perturbation theory one obtains the following 
results for either the total transition probability pit) = (1 — V(t)) or the energy spreading 
{5E{t)) 2 : 



e 2 x 



f°° duj 


i 




.(M 2 . 



F t {u)C(io) (3) 



where the expression in the square brackets should be excluded in case of SE(t). The spectral 
content of the driving is F t (w) = 2(1 - cos(wi)) for t < T/2, and F t (u) = |1 - 2c iwT / 2 + e iujt \ 
otherwise. In Fig. lb we display the rescaled energy spreading SE(t)/e for a representative 
period T and various e. For small e we observe the scaling which is implied by Eq.Q. But 
for larger e we observe a different (non-perturbative) behavior. Similar results (Fig. lc) are 
found in case of p(t)/e 2 . 

Before we go on in explaining the new "physics" in the non-perturbative regime, we would 
like to give, for sake of completeness, a very concise summary of the "physics" in the pertur- 
bative regime. The following is based on further analysis of Eq.l|3|l. extending the approach 
of Ref. [Sj . In the " standard perturbative" (or " adiabatic" ) regime, the spreading probability 
remains concentrated all the time in the initial level. In the "Wigner Lorentzian regime" the 
energy spreading profile is characterized by two energy scales and F = 2ir(ae) 2 / A. If we 
are deep in this regime we have 1 < T « Aj. The energy scale T is important for V(t), 
because this is the energy range where most of the probability is concentrated (the "core" 
component). In contrast to that the energy scale that is important for SE(t) is Af,, because 
the variance is determined by the perturbative tails of the energy distribution, and not by the 
width of the non-perturbative "core" component. This means that in the Wigner Lorentzian 
regime the "physics" of reversibility is different, depending whether we look on V(t) or on 
5E(t). To be specific, for 5E(t) we have reversibility provided t < t c i, while for V(t) we have 
reversibility for t < h/T. 

We turn now to discuss the dynamics in the non-perturbative regime, which is our main 
interest. In the absence of driving reversal [8] we get diffusion (SE(t) oc y/i) for t > t prt (e), 
where 

t ptt (e) = fi/(Vbae) . (4) 

If T/2 < t prt (e), this non-perturbative diffusion does not have a chance to develop, and 
therefore we can still trust Eq. ©. So the interesting case is T/2 > t plt (e), which means 
large enough e. It is extremely important to realize that without reversing the driving, the 
presence or the absence of E cannot be detected. It is only by driving reversal that we can 
easily determine (as in Fig. la) whether the diffusion process is reversible or irreversible. 

In order to analyze irreversibility we should distinguish between two stages in the non- 
perturbative diffusion process. The first stage (i prt < t < i s d n ) is reversible, while the second 
stage (t > i s dn) is irreversible. [For much longer time scales we have recurrences or localization 
which are not the issue of this Letter 0.] The second time scale (t s dn) is "new". It did 
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not appear in our "wavepacket dynamics" study [5], because it can be detected only by 
time/driving reversal experiment. 

We can determine the time scale i sc i n as follows. The diffusion coefficient is Z?e = 
A 2 b 5 / 2 ae/h up to a numerical prefactor [Hj. The diffusion law is SE 2 (t) = D^t. The diffusion 
process is reversible as long as E does not affects the relative phases of the participating en- 
ergy levels. This means that the condition for reversibility is (5E(t) x t)/h <C 1. The latter 
inequality can be written as t < i s dn(e), where 

(h 2 \ 1/3 ( h 3 \ 1/3 

This considerations imply that in the non-perturbative regime we should have non-universal 
scaling with respect to e 1 / 3 x T. This obviously goes beyond any implications of perturbation 
theory (Eq.JSJ)). The numerical analysis (Fig. 2) confirms the above claim. 

The non-perturbative diffusion process involves only one energy scale, and therefore in 
the non-perturbative regime the scaling of SE(T) and V(T) is similar (Fig. 3). The non 
universal scaling of V(T) versus e 1 / 3 x T is very different from the observed behavior in the 
Wigner regime (not shown) where the scaling is versus e 2 x T. Also for 8E(T) we observe the 
same non-perturbative scaling by plotting 5E(T)/e 1 / 3 versus e 1 / 3 x T. The latter should be 
contrasted with the perturbative e° x T scaling which is observed in the Wigner regime (not 
shown). The scaling of the vertical axis is explained on the basis of diffusion process after t r , 
which is of the same nature as for t < T/2. 

We are now able to identify the " sudden regime" . This is the (T, e) regime which is 
defined by the condition T <C £ s dn(e)- The condition is never satisfied in the e — > oo limit. If 
we interpret e as "velocity" then this looks quite strange in first sight, because we are used 
to associate the sudden limit with having a very "fast" process. The subtlety is resolved once 
we realize that in a proper textbook formulation of the sudden limit we should take e — > oo, 
while holding the "displacement" (e x T) equal to a constant value. 

In conclusion of the analysis of non-perturbative features of QI in the PMP sense for Wigner 
model, we would like to emphasize the distinction between "reversible" and "irreversible" 
diffusion, and the associated identification of the " sudden regime" . 
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FIG.l: Energy spreading and total transition probability as a function of time. We use units such that 
A = cr = h = 1, and b = 10. Typically we take an average over more than 100 different realizations in each 
run. (a) The actual evolution (*) with E ^ is compared with the evolution that would be obtained (o) 
if we had E = 0. The simulation is done for an extremely large perturbation (e = 1000). The remarkable 
reversibility in the E = case is an indication for the accuracy of the numerical simulations, (b) The rescaled 
energy spreading &E(t)/e, for T = 0.2, and for e < 40. The clustering of curves corresponds to e < 2.5. (c) 
The same for the rescaled total transition probability p(t)/e 2 . 



TSAMPIKOS KOTTOS, AND DORON COHEN: QUANTUM IRREVERSIBILITY 



7 




FIG. 2: (a) The time t r calculated from &E(t; T), as a function of T for amplitudes e < 40. The clustering 
of curves corresponds to e < 2.5. The inset (where only data with e > 5 arc displayed) demonstrates the 
scaling against e 1 / 3 X T in the non-perturbative regime, (b) The same, but now t r is evaluated from V{t; T). 
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FIG. 3: (a) The rescaled energy spreading SE(T)/e 1/ ' 3 versus the rescaled time e 1 / 3 X T for e values in the 
non-perturbative regime (e > 5). The inset displays the same data without rescaling. (b) The same for the 
survival probability V(T). 



